Germline homozygosity and allelic imbalance of HLA-I are common in esophagogastric adenocarcinoma and impair the repertoire of immunogenic peptides.

BACKGROUND: The individual HLA-I genotype is associated with cancer, autoimmune diseases and infections. This study elucidates the role of germline homozygosity or allelic imbalance of HLA-I loci in esophago-gastric adenocarcinoma (EGA) and determines the resulting repertoires of potentially immunogenic peptides. METHODS: HLA genotypes and sequences of either (1) 10 relevant tumor-associated antigens (TAAs) or (2) patient-specific mutation-associated neoantigens (MANAs) were used to predict good-affinity binders using an in silico approach for MHC-binding (www.iedb.org). Imbalanced or lost expression of HLA-I-A/B/C alleles was analyzed by transcriptome sequencing. FluoroSpot assays and TCR sequencing were used to determine peptide-specific T-cell responses. RESULTS: We show that germline homozygosity of HLA-I genes is significantly enriched in EGA patients (n=80) compared with an HLA-matched reference cohort (n=7605). Whereas the overall mutational burden is similar, the repertoire of potentially immunogenic peptides derived from TAAs and MANAs was lower in homozygous patients. Promiscuity of peptides binding to different HLA-I molecules was low for most TAAs and MANAs and in silico modeling of the homozygous to a heterozygous HLA genotype revealed normalized peptide repertoires. Transcriptome sequencing showed imbalanced expression of HLA-I alleles in 75% of heterozygous patients. Out of these, 33% showed complete loss of heterozygosity, whereas 66% had altered expression of only one or two HLA-I molecules. In a FluoroSpot assay, we determined that peptide-specific T-cell responses against NY-ESO-1 are derived from multiple peptides, which often exclusively bind only one HLA-I allele. CONCLUSION: The high frequency of germline homozygosity in EGA patients suggests reduced cancer immunosurveillance leading to an increased cancer risk. Therapeutic targeting of allelic imbalance of HLA-I molecules should be considered in EGA.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma.

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Frequency and functional characterization of fusion genes in squamous cell carcinoma of the lung.

INTRODUCTION: In contrast to lung adenocarcinoma (LUAD), targetable genetic alterations are less frequently detected in squamous cell carcinoma of the lung (LUSC). Over the last years, gene fusions have become promising targets in many solid cancers. Here, we analysed a cohort of LUSC, identified recurrent fusion genes and functionally characterised these tumour genomes. METHODS: A subset of 1608 squamous cell carcinomas of the lung was analysed by means of the FusionPlex® Lung Panel to identify potentially targetable gene fusions using targeted next-generation sequencing. Cases harbouring recurrent gene fusions were further analysed using FISH, Cytoscan HD arrays and cell culture experiments. RESULTS: We found both, known and novel gene fusions in about 3 % of the cases. Known fusions occurring in lung cancer included ALK::EML4, EGFRvIII, EZR::ROS1 and FGFR3::TACC. We further identified recurrent gene fusions of currently unknown biological function, involving EGFR::VSTM2A and NSD3::FGFR1 and showed that the occurrence of the EGFR::VSTM2A fusion is accompanied by high-level amplification of EGFR. Our analyses further revealed that the genomes of these LUSC patients are chromosomally unstable, which leads us to believe that such non-actionable genomic rearrangements may be a result of "chromosomal chaos" most probably not representing exclusive cancer-driving genes in this cancer entity. CONCLUSIONS: We emphasise that caution should be taken when novel fusions are found and that the appearance of new gene fusions should always be interpreted in the molecular context of the respective disease.

AI-based multi-PRS models outperform classical single-PRS models.

Polygenic risk scores (PRS) calculate the risk for a specific disease based on the weighted sum of associated alleles from different genetic loci in the germline estimated by regression models. Recent advances in genetics made it possible to create polygenic predictors of complex human traits, including risks for many important complex diseases, such as cancer, diabetes, or cardiovascular diseases, typically influenced by many genetic variants, each of which has a negligible effect on overall risk. In the current study, we analyzed whether adding additional PRS from other diseases to the prediction models and replacing the regressions with machine learning models can improve overall predictive performance. Results showed that multi-PRS models outperform single-PRS models significantly on different diseases. Moreover, replacing regression models with machine learning models, i.e., deep learning, can also improve overall accuracy.

KEAP1/NFE2L2 Pathway Signature Outperforms KEAP1/NFE2L2 Mutation Status and Reveals Alternative Pathway-Activating Mutations in NSCLC.

INTRODUCTION: Activation of the antioxidant KEAP1/NFE2L2 (NRF2) pathway leads to increased glutamine dependence and an aggressive phenotype in NSCLC. Because this pathway has been explored as a clinical target, we developed a transcriptomic signature for identifying KEAP1/NFE2L2-activated tumors. METHODS: A total of 971 NSCLC samples were used to train an expression signature (K1N2-score) to predict KEAP1/NFE2L2 mutations. There were 348 in-house NSCLCs that were analyzed using a NanoString expression panel for validation. RESULTS: The 46-gene K1N2 score robustly predicted KEAP1/NFE2L2 mutations in the validation set irrespective of histology and mutation (area under the curve: 89.5, sensitivity: 90.2%), suggesting that approximately 90% of KEAP1/NFE2L2 mutations are pathway-activating. The K1N2-score outperformed KEAP1/NFE2L2 mutational status when predicting patient survival (score p = 0.047; mutation p = 0.215). In K1N2 score-positive but KEAP1/NFE2L2 wild-type samples, enrichment testing identified SMARCA4/BRG1 and CUL3 mutations as mimics of KEAP1/NFE2L2 mutations. CONCLUSIONS: The K1N2-score identified KEAP1/NFE2L2-activated NSCLC by robustly detecting KEAP1/NFE2L2mut cases and discovering alternative genomic activators. It is a potential means for selecting patients with a constitutively active KEAP1/NFE2L2 pathway.

Dissecting the genetic heterogeneity of gastric cancer.

BACKGROUND: Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture. METHODS: We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO. FINDINGS: Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level. INTERPRETATION: Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO. FUNDING: German Research Foundation (DFG).

Why loss of Y? A pan-cancer genome analysis of tumors with loss of Y chromosome.

Loss of the Y chromosome (LoY) is frequently observed in somatic cells of elderly men. However, LoY is highly increased in tumor tissue and correlates with an overall worse prognosis. The underlying causes and downstream effects of LoY are widely unknown. Therefore, we analyzed genomic and transcriptomic data of 13 cancer types (2375 patients) and classified tumors of male patients according to loss or retain of the Y chromosome (LoY or RoY, average LoY fraction: 0.46). The frequencies of LoY ranged from almost absence (glioblastoma, glioma, thyroid carcinoma) to 77% (kidney renal papillary cell carcinoma). Genomic instability, aneuploidy, and mutation burden were enriched in LoY tumors. In addition, we found more frequently in LoY tumors the gate keeping tumor suppressor gene TP53 mutated in three cancer types (colon adenocarcinoma, head and neck squamous carcinoma, lung adenocarcinoma) and oncogenes MET, CDK6, KRAS, and EGFR amplified in multiple cancer types. On the transcriptomic level, we observed MMP13, known to be involved in invasion, to be up-regulated in LoY of three adenocarcinomas and down-regulation of the tumor suppressor gene GPC5 in LoY of three cancer types. Furthermore, we found enrichment of a smoking-related mutation signature in LoY tumors of head and neck and lung cancer. Strikingly, we observed a correlation between cancer type-specific sex bias in incidence rates and frequencies of LoY, in line with the hypothesis that LoY increases cancer risk in males. Overall, LoY is a frequent phenomenon in cancer that is enriched in genomically unstable tumors. It correlates with genomic features beyond the Y chromosome and might contribute to higher incidence rates in males.

GWAS meta-analysis of 16 790 patients with Barrett's oesophagus and oesophageal adenocarcinoma identifies 16 novel genetic risk loci and provides insights into disease aetiology beyond the single marker level.

OBJECTIVE: Oesophageal cancer (EC) is the sixth leading cause of cancer-related deaths. Oesophageal adenocarcinoma (EA), with Barrett's oesophagus (BE) as a precursor lesion, is the most prevalent EC subtype in the Western world. This study aims to contribute to better understand the genetic causes of BE/EA by leveraging genome wide association studies (GWAS), genetic correlation analyses and polygenic risk modelling. DESIGN: We combined data from previous GWAS with new cohorts, increasing the sample size to 16 790 BE/EA cases and 32 476 controls. We also carried out a transcriptome wide association study (TWAS) using expression data from disease-relevant tissues to identify BE/EA candidate genes. To investigate the relationship with reported BE/EA risk factors, a linkage disequilibrium score regression (LDSR) analysis was performed. BE/EA risk models were developed combining clinical/lifestyle risk factors with polygenic risk scores (PRS) derived from the GWAS meta-analysis. RESULTS: The GWAS meta-analysis identified 27 BE and/or EA risk loci, 11 of which were novel. The TWAS identified promising BE/EA candidate genes at seven GWAS loci and at five additional risk loci. The LDSR analysis led to the identification of novel genetic correlations and pointed to differences in BE and EA aetiology. Gastro-oesophageal reflux disease appeared to contribute stronger to the metaplastic BE transformation than to EA development. Finally, combining PRS with BE/EA risk factors improved the performance of the risk models. CONCLUSION: Our findings provide further insights into BE/EA aetiology and its relationship to risk factors. The results lay the foundation for future follow-up studies to identify underlying disease mechanisms and improving risk prediction.

Trophoblast Cell Surface Antigen 2 (TROP2) as a Predictive Bio-Marker for the Therapeutic Efficacy of Sacituzumab Govitecan in Adenocarcinoma of the Esophagus.

INTRODUCTION: The Trophoblast cell surface antigen 2 (TROP2) is expressed in many carcinomas and may represent a target for treatment. Sacituzumab govitecan (SG) is a TROP2-directed antibody-drug conjugate (ADC). Nearly nothing is known about the biological effectiveness of SG in esophageal adenocarcinoma (EAC). MATERIAL AND METHODS: We determined the TROP2 expression in nearly 600 human EAC. In addition, we used the EAC cell lines (ESO-26, OACM5.1C, and FLO-1) and a xenograft mouse model to investigate this relationship. RESULTS: Of 598 human EACs analyzed, 88% showed varying degrees of TROP2 positivity. High TROP2 positive ESO-26 and low TROP2 positive OACM5.1C showed high sensitivity to SG in contrast to negative FLO-1. In vivo, the ESO-26 tumor shows a significantly better response to SG than the TROP2-negative FLO-1 tumor. ESO-26 vital tumor cells show similar TROP2 expression on all carcinoma cells as before therapy initiation, FLO-1 is persistently negative. DISCUSSION: Our data suggest that sacituzumab govitecan is a new therapy option in esophageal adenocarcinoma and the TROP2 expression in irinotecan-naïve EAC correlates with the extent of treatment response by sacituzumab govitecan. TROP2 is emerging as a predictive biomarker in completely TROP2-negative tumors. This should be considered in future clinical trials.

Genome instability is associated with ethnic differences between Asians and Europeans in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the deadliest cancer types with diverse etiological factors across the world. Although large scale genomic studies have been conducted in different countries, integrative analysis of HCC genomes and ethnic comparison across cohorts are lacking. Methods: We first integrated genomes of 1,349 HCC patients from five large cohorts across the world and applied multiple statistical methods in identifying driver genes. Subsequently, we systematically compared HCC genomes and transcriptomes between Asians and Europeans using the TCGA cohort. Results: We identified 29 novel candidate driver genes, many of which are infrequent tumor suppressors driving late-stage tumor progression. When we systematically compared ethnic differences in the genomic landscape between Asian and European HCCs using the TCGA cohort (n = 348), we found little differences in driver frequencies. Through multi-modal integrative analysis, we found higher genomic instability in Asians together with a collection of molecular events ranging from tumor mutation burden (TMB), copy number alterations as well as transcriptomic subtypes segregating distinctively between two ethnic backgrounds. Strikingly, we identified an Asian specific transcriptomic subtype with multiple ethnically enriched genomic alterations, in particular chromosome 16 deletion, leading to a clinically aggressive RNA subgroup unique to Asians. Integrating multi-modal information, we found that survival models predict patient prognosis much better in Asians than in Europeans, demonstrating a higher potential for precision medicine applications in Asia. Conclusion: For the first time, we have uncovered an unprecedented amount of genomic differences segregating distinctively across ethnicities in HCC and highlighted the importance of differential disease biology and management in HCC across ethnic backgrounds.

Observations that suggest a contribution of altered dermal papilla mitochondrial function to androgenetic alopecia.

Androgenetic alopecia (AGA) is a prevalent hair loss condition in males that develops due to the influence of androgens and genetic predisposition. With the aim of elucidating genes involved in AGA pathogenesis, we modelled AGA with three-dimensional culture of keratinocyte-surrounded dermal papilla (DP) cells. We co-cultured immortalised balding and non-balding human DP cells (DPCs) derived from male AGA patients with epidermal keratinocyte (NHEK) using multi-interfacial polyelectrolyte complexation technique. We observed up-regulated mitochondria-related gene expression in balding compared with non-balding DP aggregates which indicated altered mitochondria metabolism. Further observation of significantly reduced electron transport chain complex activity (complexes I, IV and V), ATP levels and ability to uptake metabolites for ATP generation demonstrated compromised mitochondria function in balding DPC. Balding DP was also found to be under significantly higher oxidative stress than non-balding DP. Our experiments suggest that application of antioxidants lowers oxidative stress levels and improves metabolite uptake in balding DPC. We postulate that the observed up-regulation of mitochondria-related genes in balding DP aggregates resulted from an over-compensatory effort to rescue decreased mitochondrial function in balding DP through the attempted production of new functional mitochondria. In all, our three-dimensional co-culturing revealed mitochondrial dysfunction in balding DPC, suggesting a metabolic component in the aetiology of AGA.

Immune responses against shared antigens are common in esophago-gastric cancer and can be enhanced using CD40-activated B cells.

BACKGROUND: Specific immune response is a hallmark of cancer immunotherapy and shared tumor-associated antigens (TAAs) are important targets. Recent advances using combined cellular therapy against multiple TAAs renewed the interest in this class of antigens. Our study aims to determine the role of TAAs in esophago-gastric adenocarcinoma (EGA). METHODS: RNA expression was assessed by NanoString in tumor samples of 41 treatment-naïve EGA patients. Endogenous T cell and antibody responses against the 10 most relevant TAAs were determined by FluoroSpot and protein-bound bead assays. Digital image analysis was used to evaluate the correlation of TAAs and T-cell abundance. T-cell receptor sequencing, in vitro expansion with autologous CD40-activated B cells (CD40Bs) and in vitro cytotoxicity assays were applied to determine specific expansion, clonality and cytotoxic activity of expanded T cells. RESULTS: 68.3% of patients expressed ≥5 TAAs simultaneously with coregulated clusters, which were similar to data from The Cancer Genome Atlas (n=505). Endogenous cellular or humoral responses against ≥1 TAA were detectable in 75.0% and 53.7% of patients, respectively. We found a correlation of T-cell abundance and the expression of TAAs and genes related to antigen presentation. TAA-specific T-cell responses were polyclonal, could be induced or enhanced using autologous CD40Bs and were cytotoxic in vitro. Despite the frequent expression of TAAs co-occurrence with immune responses was rare. CONCLUSIONS: We identified the most relevant TAAs in EGA for monitoring of clinical trials and as therapeutic targets. Antigen-escape rather than missing immune response should be considered as mechanism underlying immunotherapy resistance of EGA.

Variant profiling of colorectal adenomas from three patients of two families with MSH3-related adenomatous polyposis.

The spectrum of somatic genetic variation in colorectal adenomas caused by biallelic pathogenic germline variants in the MSH3 gene, was comprehensively analysed to characterise mutational signatures and identify potential driver genes and pathways of MSH3-related tumourigenesis. Three patients from two families with MSH3-associated polyposis were included. Whole exome sequencing of nine adenomas and matched normal tissue was performed. The amount of somatic variants in the MSH3-deficient adenomas and the pattern of single nucleotide variants (SNVs) was similar to sporadic adenomas, whereas the fraction of small insertions/deletions (indels) (21-42% of all small variants) was significantly higher. Interestingly, pathogenic somatic APC variants were found in all but one adenoma. The vast majority (12/13) of these were di-, tetra-, or penta-base pair (bp) deletions. The fraction of APC indels was significantly higher than that reported in patients with familial adenomatous polyposis (FAP) (p < 0.01) or in sporadic adenomas (p < 0.0001). In MSH3-deficient adenomas, the occurrence of APC indels in a repetitive sequence context was significantly higher than in FAP patients (p < 0.01). In addition, the MSH3-deficient adenomas harboured one to five (recurrent) somatic variants in 13 established or candidate driver genes for early colorectal carcinogenesis, including ACVR2A and ARID genes. Our data suggest that MSH3-related colorectal carcinogenesis seems to follow the classical APC-driven pathway. In line with the specific function of MSH3 in the mismatch repair (MMR) system, we identified a characteristic APC mutational pattern in MSH3-deficient adenomas, and confirmed further driver genes for colorectal tumourigenesis.

Complementary Sequential Circulating Tumor Cell (CTC) and Cell-Free Tumor DNA (ctDNA) Profiling Reveals Metastatic Heterogeneity and Genomic Changes in Lung Cancer and Breast Cancer.

INTRODUCTION: Circulating tumor cells (CTCs) and cell-free tumor DNA (ctDNA) are tumor components present in circulation. Due to the limited access to both CTC enrichment platforms and ctDNA sequencing in most laboratories, they are rarely analyzed together. METHODS: Concurrent isolation of ctDNA and single CTCs were isolated from lung cancer and breast cancer patients using the combination of size-based and CD45-negative selection method via DropCell platform. We performed targeted amplicon sequencing to evaluate the genomic heterogeneity of CTCs and ctDNA in lung cancer and breast cancer patients. RESULTS: Higher degrees of genomic heterogeneity were observed in CTCs as compared to ctDNA. Several shared alterations present in CTCs and ctDNA were undetected in the primary tumor, highlighting the intra-tumoral heterogeneity of tumor components that were shed into systemic circulation. Accordingly, CTCs and ctDNA displayed higher degree of concordance with the metastatic tumor than the primary tumor. The alterations detected in circulation correlated with worse survival outcome for both lung and breast cancer patients emphasizing the impact of the metastatic phenotype. Notably, evolving genetic signatures were detected in the CTCs and ctDNA samples during the course of treatment and disease progression. CONCLUSIONS: A standardized sample processing and data analysis workflow for concurrent analysis of CTCs and ctDNA successfully dissected the heterogeneity of metastatic tumor in circulation as well as the progressive genomic changes that may potentially guide the selection of appropriate therapy against evolving tumor clonality.

Integrative Profiling of T790M-Negative EGFR-Mutated NSCLC Reveals Pervasive Lineage Transition and Therapeutic Opportunities.

PURPOSE: Despite the established role of EGFR tyrosine kinase inhibitors (TKIs) in EGFR-mutated NSCLC, drug resistance inevitably ensues, with a paucity of treatment options especially in EGFR T790M-negative resistance. EXPERIMENTAL DESIGN: We performed whole-exome and transcriptome analysis of 59 patients with first- and second-generation EGFR TKI-resistant metastatic EGFR-mutated NSCLC to characterize and compare molecular alterations mediating resistance in T790M-positive (T790M+) and -negative (T790M-) disease. RESULTS: Transcriptomic analysis revealed ubiquitous loss of adenocarcinoma lineage gene expression in T790M- tumors, orthogonally validated using multiplex IHC. There was enrichment of genomic features such as TP53 alterations, 3q chromosomal amplifications, whole-genome doubling and nonaging mutational signatures in T790M- tumors. Almost half of resistant tumors were further classified as immunehot, with clinical outcomes conditional on immune cell-infiltration state and T790M status. Finally, using a Bayesian statistical approach, we explored how T790M- and T790M+ disease might be predicted using comprehensive genomic and transcriptomic profiles of treatment-naïve patients. CONCLUSIONS: Our results illustrate the interplay between genetic alterations, cell lineage plasticity, and immune microenvironment in shaping divergent TKI resistance and outcome trajectories in EGFR-mutated NSCLC. Genomic and transcriptomic profiling may facilitate the design of bespoke therapeutic approaches tailored to a tumor's adaptive potential.